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Abstract 

We study big bang nucleosynthesis in the presence of large mass-scale, non-linear 
entropy fluctuations. Overdense regions, with masses above the local baryon- Jeans mass, 
are expected to collapse and form condensed objects. Surviving nucleosynthesis products 
therefore tend to originate from underdense regions. We compute expected surviving 
light element (^H, '^He, ^He, ^Li) abundance yields for a variety of stochastic fluctuation 
spectra. In general, we find that spectra with significant power in fiuctuations on length 
scales below that of the local baryon Jeans mass produce nucleosynthesis yields which are 
in confiict with observationally inferred primordial abundances. However, when this small 
scale structure is absent or suppressed, and the collapse efficiency of overdense regions is 
high, there exists a range of fiuctuation spectral characteristics which meet all primordial 
abundance constraints. In such models abundance constraints can be met even when 
the pre-coUapse baryonic fraction of the closure density is Qb ~ 0.2h~^ {h is the Hubble 
parameter in units of 100 km s~^Mpc~^). Nucleosynthesis in these models is characterized 
by high ^H/H and low ^He mass fraction relative to a homogeneous big bang at a given 
value of Q^^h?. A potentially observable signature of these models is the production of 
intrinsic primordial abundance variations on baryon mass-scales up to IQ^^Mq — lO^^M©. 

Subject headings: cosmology: theory - early universe - 
theory - nucleosynthesis, abundances - theory - large scale structure: general 
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1. Introduction 

In this paper we calculate the nucleosynthesis to be expected in the presence of primor- 
dial isocurvature baryon number (hereafter, PIB) fluctuations. Such fluctuations recently 
have been proposed as possible seeds for large scale structure formation (Peebles 1987ab; 
Suginohara & Suto 1992; Cen, Ostriker, & Peebles 1993). Density fluctuations in PIB mod- 
els are essentially entropy fluctuations. In an earlier series of papers (Jedamzik & Fuller 
1994; Jedamzik, Fuller, & Mathews 1994) we have examined the evolution and nucleosyn- 
thesis effects of small mass-scale nonlinear entropy fluctuations. In the present paper we 
extend our nucleosynthesis study to larger mass-scale fluctuations, including those relevant 
for PIB models. 

Proposed PIB models invoke a spectrum of entropy fluctuations characterized by in- 
creasing amplitude with decreasing fluctuation mass scale. In these models the relationship 
between the fluctuation amplitude, 5p/ p , and the mass contained inside the region of the 
fluctuation, M , can be written as 

^~M-^-t, (1) 
P 

where n is a spectral index. Spectra of this form are fully specified by giving the index, 
n , and the mass scale, Mu^ity ; on which Sp/p = 1. Typical spectral indices employed in 
PIB models include n = Oorn=— 1, where accompanying values of Munity are 

IQii > > io6 . (2) 

We have no compelling microscopic theory for how PIB fluctuations could be generated 
in general, though there are some cogent suggestions. Yokoyama & Suto (1991) and Dolgov 
& Silk (1992) have proposed microscopic mechanism for the generation of baryon number 
(entropy) fluctuations. These schemes exploit quantum fluctuations in a CP-violating 
angle to effect a spatial variation in baryon-to-photon number. The resulting fluctuations 
can be on the large mass-scales of interest in PIB models if the baryon-generating process 
occurs before, or during, an inflationary epoch. 

A common, though not necessarily inevitable, feature of inflationary baryogenesis 
models is the production of both matter and antimatter domains separated by astrophysi- 
cally signiflcant length scales. The production of antimatter domains can be avoided if the 
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baryon number fluctuation amplitudes are below unity on all length scales. In this case, 
however, the resulting light element abundance yields from nucleosynthesis will be only 
slightly altered from abundance yields in homogeneous big bang nucleosynthesis at the 
same average baryon-to-photon ratio rj (for a discussion of small amplitude fluctuations in 
r] and primordial nucleosynthesis yields see Epstein & Petrosian 1973). 

Though large scale structure considerations in PIB models utilize only the large mass- 
scale linear {6p/p < 1) end of the spectrum in equation (1), successful models may well 
require some small mass-scale, nonlinear {Sp/p > 1) structure in order to produce coUi- 
sionless dark matter and/or to effect late re-ionization of the universe (Cen et al. 1993). 
Of course, the biggest hurdle which PIB models face may well be stringent observational 
limits on cosmic background radiation anisotropy (Gorski & Silk 1989; Chiba, Sugiyama, 
& Suto 1993; Hu & Sugiyama 1994). 

If any PIB model for structure formation must invoke a nonlinear lower mass-scale 
fluctuation tail in a spectrum like that in equation (1), or if such a tail is an inevitable 
consequence of some microscopic fluctuation generation mechanism, then primordial nu- 
cleosynthesis effects may well provide constraints on the model. These constraints could be 
complimentary to cosmic-background-anisotropy constraints, as nucleosynthesis probes a 
different end of the fluctuation spectrum. In this paper we aim to see under what circum- 
stances such constraints on PIB- like models could be found. Our study and our derived 
constraints will generally apply to any models which have nonlinear entropy fluctuations - 
not just PIB-likc models for large scale structure formation. 

Studies of homogeneous big bang nucleosynthesis (HBBN) only allow a narrow range 
for the fraction of the closure density that can be contributed by baryons 

0.01 < VL^h^ < 0.015 , (3) 

where h is the Hubble parameter in units of 100 kms~^Mpc~^ (cf. Walker et al. 1991; 
Smith, Kawano, & Malaney 1993). Models of inhomogeneous big bang nucleosynthesis 
(IBBN) with sub-horizon scale nonlinear entropy fluctuations cannot significantly change 
the conlusions of HBBN as to the upper limit on ri^/i^ (cf. Jedamzik, Fuller, & Mathews 
1994; Thomas et al. 1994). 

However, there has been a longstanding suggestion that large mass scale entropy fluc- 
tuations might provide a way to circumvent the HBBN upper bound on Q.^^h'^ (Harrison 
1968; Zeldovich 1975; Rees 1984; Sale & Mathews 1986). In these schemes regions of 
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relatively large baryon-to-photon ratio (rj) collapse to form inert remnants and therefore 
remove their "bad" nucleosynthesis products from the primordial medium. The nucleosyn- 
thesis products which survive in these models reflect the relatively lower baryon density 
of regions which do not collapse. In principle, the pre-coUapse flbh^ could be considerably 
larger than the value of this quantity (presumbly O^/i^ ^ 0.01 — 0.015) which characterizes 
the low density medium which surrounds collapsed regions. 

Indeed, studies of the fate of large mass entropy fluctuations in the epoch between 
primordial nucleosynthesis and some time after recombination indicate that overdense re- 
gions would collapse with high efficiency (Hogan 1978, 1993; Kashlinsky & Rees 1983; Loeb 
1993). These studies show that the evolution of large entropy fluctuations is dominated by 
photon-electron Thomson drag (Hogan 1993; Loeb 1993). Any fluctuation with a baryon 
mass larger than about 



will shrink slowly, with Thomson drag providing an efficient mechanism for cooling the col- 
lapsing protons and electrons, and also damping rotation. The end result of this process 
will be the production of a condensed object, either a black hole or possibly, if fragmen- 
tation occurs, many small brown dwarfs (Kashlinsky & Rees 1983; Hogan 1993). Note 
that the mass scale in equation (4) is essentially coincident with the local baryon Jeans 
mass. Of course, the true Jeans mass in a radiation dominated environment is close to the 
horizon mass. By local baryon Jeans mass we mean the effective Jeans mass one would 
calculate inside a fluctuation by neglecting all photon stresses. 

We note, however, that the studies of Harrison (1968), Zeldovich (1975), Rees (1984), 
and Sale &; Mathews (1986) provide what are at best simplistic treatments in their calcu- 
lations of nucleosynthesis yields in the presence of collapsing regions. In fact, all of these 
studies except Sale & Mathews (1986) treat the distribution of 77 in a two- phase, bi-modal 
fashion. All of these studies neglect the significant complication that fiuctuations may 
reside inside larger mass-scale fiuctuations, the so called cloud-in-cloud problem. None of 
these prior studies addresses the nucleosynthesis effects of fiuctuations below that in equa- 
tion (4) - such fluctuations are damped by expansion against Thomson drag (cf. Hogan 
1993; Jedamzik & Fuller 1994), but on time scales that allow signiflcant nucleosynthesis 
effects (Alcock et al. 1990; Jedamzik, Fuller, & Mathews 1994). 




(4) 
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Recently, Gnedin, Ostriker, & Rees (1994) have re-examined the problem of nucleosyn- 
thesis with collapsing entropy fluctuations. These authors provide a more sophisticated 
numerical treatment of nucleosynthesis from regions of varying density and, additionally, 
attempt to take account of light element "reprocessing" effected by accretion on black holes 
(Gnedin &: Ostriker 1992). They conclude that significant power on fiuctuation scales be- 
low the limit in equation (4) could lead to unacceptable nucleosynthesis which, in turn, 
could wreck any attempt to employ collapsing high-mass-scale fluctuations to circumvent 
HBBN bouns on fl^h'^. They do, however provide for several loopholes in this conclusion: 
phase correlations in fluctuations; and invocation of flnely tuned IBBN scenarios for small 
mass-scale fluctuations. 

However, Gnedin, Ostriker, & Rees (1994) do not attempt to take account of the 
effects of the cloud-in-cloud problem. Furthermore, their suggestion that baryon diffusive 
effects on small scales may provide a loophole on the il^h^ limit is suspect, as it requires 
an extreme, though not impossible, degree of fine tuning. 

In this paper we attempt a detailed numerical treatment of nucleosynthesis with vari- 
ous stochastic entropy fluctuation spectra. We explicitly treat the cloud-in-cloud problem. 
We flnd that any signiflcant power in stochastic fluctuations on mass-scales smaller than 
that in equation (4) inevitably leads to overproduction of ^He and ^Li relative to obser- 
vationally inferred limits on these abundances. Furthermore, we argue that persistence of 
fluctuations down to mass scales on which baryon diffusive effects are signiflcant is unlikely 
to change our conclusions. Finally, we show that only under a restrictive set of circum- 
stances (i.e. no signiflcant power in fluctuations on mass scales below Mj and high collapse 
efliciency for scales above Mj) it is possible to evade HBBN bounds on Vt^h? . We note 
that evasion of the bound on Q.i,h'^ in these models also usually requires that the primordial 
abundance of ^Li be close to the Population I vaule of ^Li/H ^ 10~^. 

In Section 2 we discuss our model for fluctuation evolution and nucleosynthesis and 
the numerical techniques employed. In Section 3 we present results for nucleosynthesis 
yields for various fluctuation spectra and assumed fluctuation evolution parameters. We 
compare these to observationally inferred abundances. We give conclusions in Section 4. 
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2. Simulations of baryon-to-photon fluctuations 

In this section we discuss our numerical calculations of fluctuation evolution and 
nucleosynthesis. The basis of our models is a stochastic distribution of fluctuations in 
baryon-to-photon ratio on various scales. These scales will be selected in such a manner as 
to roughly approximate the density distribution in an overall spectrum like that in equation 
(1). Note, however, that much is hidden in a simplistic power law density distribution 
like equation (1). Going from equation (1) to a numerical representation of fluctuation 
amplitude and mass-scale distribution is not a unique procedure. In what follows we 
present our numerical approach to this problem and discuss the statistics of our numerically 
generated distributions of baryon-to-photon number. 

2.1. The numerical models 

As a first step, we generate a stochastic distribution of baryon-to-photon number 
by employing a spatially inhomogeneous gaussian random variable, A{x). This gaussian 
random variable can be Fourier decomposed to give, 

A(x) = 2^ AkCos(kx + (f)k) , (5) 
fc 

where a; is a spatial coordinate. In this expression the amplitudes A^ and phases (f)k are 
chosen randomly according to the distribution functions, 

PMk) = ^ 0fee[O,27r], (6a) 

1 1 f 1 A^\ 
PA(^fe) = exp --^ ylfe e [-00, oo] . (66) 

V27r CTfc \ 2 a^J 

Note that our adoption of the distribution function in equation (6a) implies that there 
will be no correlations between the phases of different Fourier modes - of course, this 
is the definition of a gaussian random variable. We adopt a power law dependence on 
wavenumber for the variance, ak, 

. (7) 
where n is a spectral index, not to be confused with that in equation (1). 
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The spatially dependent baryon-to-photon ratio , r]{x) , is taken to be a function of 
the spatially fluctuating random variable A{x). We consider three different models: 



ri{x) = r]NA'^{x) ; 



(8a) 



ri{x) = riNA^^{x) ; 



{8b) 



r]{x) =?7iv(A^°(x) +a) ; 



(8c) 



where tjn and a are constants. Note that the functions in equations (8abc) were deliberately 
chosen to be positive-definite. This choice guarantees that antimatter-domains are avoided. 

Our choice of the spatial distributions in equations (8abc) is not based on specific 
baryogenesis scenarios. Rather, we have chosen these models simply to generate a wide 
variety of stochastic baryon number distributions. It remains to be shown if baryon number 
distributions of these, or similiar characteristics, could arise naturally during the evolution 
of the very early universe. 

In this paper we analyze a one-dimensional analogue to a three-dimensional distri- 
bution of entropy (or, equivalently, rj). By using a one-dimensional distribution of rj we 
are able to investigate a much wider range of mass scales in a single simulation. The 
analogy to the three-dimensional theory is attained by replacing the spectral index n in 
equation (7) with n/3. For example, n = —3 corresponds to a scale- invariant Harrison- 
Zeldovich spectrum in the three-dimensional theory; whereas n = —1 corresponds to the 
Harrison- Zeldovich spectrum in one dimension. 

In what follows we will always refer to the three-dimensional spectral index when we 
wish to characterize a particular one-dimensional distribution in rj. In a one-dimensional 
calculation we compute the "mass" of a region by multiplying "density" (or amplitude 
5p/p) by the length scale of the region. In three dimensions, masses are the product of 
Sp/p and a volume. The transformation of spectral indices n ^ ?^/3, insures that the 
functional dependence of 6p/p on mass-scale (e.g. equation 1) that we derive in our one- 
dimensional calculation is the same as that in an equivalent (n/3 >—>■ n) three-dimensional 
case. 

It is well known that diffusive and hydrodynamic processes occuring during the epoch 
of primordial nucleosynthesis can significantly alter the light-element abundance yields. 
Diffusive processes will play an important role during nucleosynthesis when significant 
small-scale fluctuations on mass scales between M ~ lO^^-'^M© and M ~ 1{)~^^Mq are 
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present. Nuclear abundance yields in such scenarios will certainly be highly geometry 
dependent so that a full three-dimensional treatment would be needed. (Mathews et al. 
1990; Meyer et al. 1991). However, we will suggest below that, for a truly stochastic baryon 
number distribution (e.g. minimal phase correlations between different Fourrier modes of 
the distribution) , the presence of such small-scale fluctuations cannot be reconciled with the 
observationally inferred primordial abundance limits. We will therefore not be concerned 
with diffusive processes. 

In the case where diffusive and hydrodynamic processes during nucleosynthesis are 
unimportant, average light-clement abundance yields arc determined by a weighted average 
over the standard homogeneous big bang yields of separate regions at different baryon-to- 
photon ratios. In this study we obtain the light-element abundance yields by employing 
the homogeneous big bang nucleosynthesis code of Wagoner, Fowler, & Hoyle (1967) as 
updated by Kawano (1992). 

In our calculations the light-element contributions from overdense regions, which are 
expected to collapse and form compact objects, are excluded from the abundance average. 
We choose all regions which are overdense on average by some critical amount Acr to be 
"destined" for collapse. We define the overdensity parameter. Ax , to be. 



In the above expression < -q > denotes the cosmic average baryon-to-photon ratio and 
r]x (x) represents the average baryon-to-photon ratio within a region of size A around space 
coordinate x. In our numerical prescription only regions with = Acr are taken to 
collapse. If a region is found to have A;^ > A^r we look for a larger region, which contains 
the original, for which A> = Acr- This larger region is taken to collapse. In other words, 
we take the largest scale for which Ax — A^r to collapse and take everything inside to 
oblivion. 

The above prescription for determining which regions collapse is a widely used analytic 
tool in the study of structure formation (cf. Press & Schechter 1974). This procedure is 
also confirmed by comparison of the analytic results to numerical simulations (Efstathiou 
et al. 1988). Note that the Press-Schechter analysis has been recently modified in order 
to account properly for the cloud- in-cloud problem (Jedamzik 1994). 

The epoch at which an overdense region collapses is determined, for a given cosmolog- 
ical model, predominantly by its initial average overdensity. For example, the formation of 




(9) 



<r]> 
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supermassive black holes in PIB models requires early collapse and high initial overdensities 
(Hogan 1993). 

Partly, this requirement arises because Thomson-drag on background photons be- 
comes a less efficient mechanism for entropy transport and rotational damping at lower 
cosmic temperature. A certain determination of the fate of a collapsing cloud would re- 
quire a full three-dimensional hydrodynamic simulation. Determining the fate of overdense 
regions is important, since we need to know whether or not to count their freeze-out big 
bang nucleosynthesis products in the final abundance yields of the diffuse, "uncoUapsed" 
, surviving background of primordial gas. 

The value of Agj. where the efficiency of collapse becomes large depends on the cos- 
mological model and the ionization history of the universe (Hogan 1993). We define the 
collapse efficiency parameter, / , to be the fraction of regions with A a = Acr which actu- 
ally do collapse. In this paper we will simply treat Acr and / as parameters to be varied. 
We will determine the sensitivity of our calculated light-element abundance yields, and 
ratio of diffuse baryons to dark baryons, to these parameters. 

To determine the light-element abundances in the presence of fluctuations it is not 
sufficient to consider only the average densities of regions. Rather, in principle, a detailed 
knowledge of the baryon number distribution on all scales is required. Toward meeting 
this requirement, we define Pr){x) to be the probability that the region at point x has 
baryon-to- photon ratio rj. We will find, ultimately, that all small mass scale fiuctuations 
(M < Mj) would have to be suppressed in order that a PIB-like fluctuation spectrum 
would be able to produce nucleosynthesis consistent with observational constraints. If we 
utilize this result, and simply assume that on small enough scales fluctuation amplitudes 
go to zero, then we can exploit the overall translational symmetry of the universe to write 
Priix) = P{r]). In other words the probability to flnd a region with baryon-to-photon ratio 
T] is independent of x; mathematically Pij]) is the probability of flnding any point to have 
baryon-to-photon ratio rj. By translational symmetry in this argument we mean that any 
one large region of the universe (large enough to have ff) must be fully equivalent and 
similar to any other such region. 

Note that writing P{ri) makes sense only if there is a small scale cutoff in structure. 
In our numerical calculations any scale smaller than this cutoff scale is taken to be homo- 
geneous. The probability for flnding any such spatial zone to have baryon-to-photon ratio 
?7 is Pirj). The probability that a given scale A has average baryon-to-photon ratio 77 is 
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defined to be P\{r]). Clearly, with this definition, Px{r]) = P{r]) for all A < Ac, where Ac 
is the cutoff scale. We will discuss the statistical relationship between P\{r]) and P{r]) for 
larger scales below. 

It is quite important for what follows to note that P{r]), and any mass-based statis- 
tical criteria for describing mass distributions, will be invariant under the transformation 
from one dimension to three dimensions as outlined above. An example of an invariant 
mass-based distribution function is the following: the distribution of masses which have 
overdensities equal to Acr- 

We wish to stress the necessity of a numerical treatment for the reliable determina- 
tion of nuclear abundance yields. Furthermore, there are many subtle pitfalls involved in 
estimating abundance yields in an inhomogeneous environment. For example, it is not ad- 
equate to simply average over nucleosynthesis yields resulting from all regions below some 
threshold, r] < rjcr ■ This procedure is inadequate since such regions could possibly be 
underdense clouds within collapsing overdense regions. As another example, an analytic 
computation of Px{r]) is difficult (impossible) since the transformation in equations (8abc) 
introduces phase correlations between the different Fourier modes of the baryon number 
distribution. We therefore expect the stochastic baryon number distribution of our model 
to have a non-gaussian character. 

We have simulated three types of stochastic baryon number distributions by applying 
the transformation in equation (8abc) to a gaussian random variable. For these simulations 
we generated 10^ Fourier modes of the gaussian random variable. The wavevectors of 
these Fourier modes were determined by applying periodic boundary conditions in a box 
extending from zero to one. A small part of these distributions are shown in Figures labc. 
Note that the value of 10° on the logarithmic abscissa corresponds to the cosmic average 
baryon-to-photon ratio, < 77 > , in all three figures. 

In Figure la most of the peaks in the distribution of 77 are seen to be between values of 
5 and 10 (e.g., five to ten times average density); whereas, the distribution shown in Figure 
lb exhibits peaks with values between 10 and 100. The distribution of i] in Figure Ic also 
exhibits large overdense peaks. However, in this distribution there are no regions with very 
low density. The dotted boxes in these figures indicate the regions which are overdense 
by the critical amount, Acr- These regions will be expected to collapse very early on. 
It is obvious from the figures that such overdense regions often are centered around very 
prominent peaks. Note, however, that these overdense regions can include underdense 
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material at small baryon-to-photon ratios as well. We will analyze the nucleosynthesis 
resulting from these baryon-to-photon number distributions in Section 3. 



2.2. Statistics of the baryon-to-photon number distributions 

In this section wc will investigate the statistics of the baryon-to-photon number dis- 
tributions generated by the prescription outlined in section 2.1. Consider first a gaussian 
random variable A. The probability distribution P{A) for finding a small region to have 
value A is given by the gaussian distribution, 

11 / 1 A 

-= expl 

V27r (Jtot V 2 cT^ot 



PiA) = ^-_exp( --^ ) , (10a) 



where 



1/ L\ 



Note the relationship between these expressions and those for P/i(Afc) and af, in equations 
(6b) and (7). Here (L/27r)~^ is the volume in wavevector space in which there is one Fourier 
mode. Equation (10) makes use of the fact that the sum, A^ of several normally distributed 
quantities (e.g., the amplitudes of the uncorrelated Fourier modes, A]^) itself follows a 
gaussian distribution. The square of the variance, a^Q^ , of this gaussian distribution in A 
is then given by the sum of the squares of the individual variances of the Fourier modes, 
cr^. Note that in equation (10b) we introduce a cutoff- wavevector kc-, since we assume that 
there are no small-scale fluctuations. In particular, we take ak ^ ^ ior k > kc = 27r/Ac. 

With the help of equation (10a) we can derive the probability function P{ri). Assuming 
a transformation of the form r]{x) — r]N{A'^'^{x) + a) we obtain, 

V 27r matotVN Kvn J { cr^ot ) 

where r] falls in the range ry e [ar/jv, oo]. In this expression note that m is an integer index, 
while r/jv and a are constants. 

It should be stressed that it is not possible to obtain the probability function P\{'q) 
by simply employing equation (11) with a modified variance atot i— > cr ~ dk'a'^, . This is 
because with a positive definite transformation such as that in equation (8abc) the average 
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rj of a. region is determined by the larger-scale Fourier modes with k' < k as well as by the 
smaller-scale Fourier modes k' > k. We therefore must resort to numerical techniques in 
order to analyze the baryon-to-photon number distribution on larger scales. 

We define the integrated variance in the distribution of r] on mass scales below the cut- 
off scale to be cr(M < Mc). If Mc is the mass scale corresponding to the cutoff wavevector 
kc, then we can define 



1 

— — / dri{ri- < rj >fP{r]) \ . 



(12) 



The quantity cr^ determines the likely magnitude of fluctuation amplitudes. For the trans- 
formation 'q{x) = rjN{-A^^{x) + aj we can derive, for example, 



V^r(2m+^)-r2(m + ^) 



(13) 



where the average baryon-to-photon ratio is given by 



<r]>-- 



(2ay™r(m+i) + ar(i) 



(14) 



If we apply equation (13) to the transformation in equation (8a) we obtain ar^{M < Mc) = 
V2. By contrast, the fluctuation amplitudes for the transformations in equations (8bc) are 
larger and yield arf{M < Mc) ~ 27. It should be noted that cr^ is independent of atot- 

We have investigated numerically the statistics of the baryon-to-photon number dis- 
tributions on large scales. Figures 2abc show the distribution functions Pxiv) on different 
mass scales M (or, equivalently, spatial scales A) for the three transformations in equa- 
tions (Babe). In these figures the value of 1 on the abscissa corresponds to the average 
baryon-to-photon ratio; whereas, the scale of the ordinate is chosen in each case so that 
the complete distribution of P\{rj) is displayed. The lines on these figures are for the dis- 
tribution function Px{r]) on mass scales M = Mc (sofid), M = lOM^ (dotted), M = lOOMc 
(short-dashed), and M — lOOOMc (long-dashed). We will determine the numerical value 
of Mc in the following section. The dashed-dotted line shows the result from equation (11) 
for comparison. It can be seen from these figures that the distribution functions Px{r]) 
are highly non-gaussian on small scales, but approach a gaussian distribution character on 
large scales. Note that the variance, , decreases as the mass scale increases. 
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A non-gaussian distribution P\{r]) on small scales could only result if there were 
phase correlations between the different Fourier modes of the baryon-to-photon number 
distribution. We note that phase correlations are introduced in our procedure as a result of 
our application of the transformation equation (8abc) to the uncorrelated gaussian random 
variable. Similiar results (e.g., non-gaussian distributions on small scales and gaussian 
distributions on large scales) have been obtained by Yamamoto et al. (1992). These 
authors analyzed the baryon-to-photon number distribution resulting from the inflationary 
baryogenesis scenario proposed by Yokoyama & Suto (1999). That baryogenesis scenario 
assumed that the production of baryon number was proportional to the a trigonometric 
function of a gaussian random variable, e.g. 'r]{x) = r]Qsin{A[x)) . 

In Figures 3abc we show the fluctuation variance a^{M) as a function of mass scale 
for the three different transformations in equations (8abc). In each of these flgures the 
dotted line gives cr^(M) in the case where the variance of the Fourier amplitudes is taken 
to be constant, cr^ = constant. The dashed line in these figures corresponds to choosing 
cr| ~ k~'^-^. Here, the index —2.4 is the "three-dimensional" index. For purposes of 
comparison in these figures we show the root-mean-square mass fluctuation expected for 
a gaussian random variable, 

_ 1 n 

The solid lines in these flgures correspond to ar^{M) from equation (15) with n = (lower 
line) and n= —1.5 (upper line). 

2.3. Normalization of the baryon-to-photon number distributions 

We can now determine the absolute magnitude of the cutoff mass-scale, , for 
models with different spectral indices , n , and variances, (7^(Mc). Primordial baryon- 
to-photon number fluctuations are constrained by the observed structure of the universe 
and by observationally inferred primordial light element abundances. Fluctuations are also 
constrained by the high degree of isotropy of the cosmic microwave background radiation 
(hereafter, CMBR). However, we will assume here that an early reionization can erase any 
preexisting anisotropics in the CMBR caused by the baryon-to-photon number fluctuations. 
It is not at all obvious whether a given PIE spectrum of fluctuations will result in early 
reionization and thus have a chance of escaping constraints from CMBR anisotropy limits. 



13 



The variance, or root-mean-square fluctuations in mass, (5M/M)r.m.s. ) at the present 
epoch has been determined by Davis & Peebles (1983). They find (5M/M)r.m.s. ~ 1 in 
a volume of size (8h~-'^Mpc)^. This volume would correspond to a mass scale of Mg 
6 X 10^^ MQ{ilhh^ /0.1)h~^ . It is well known that sub-horizon, super-Jeans-mass size baryon 
number perturbations grow proportionally to the scale factor of the universe during a 
matter dominated epoch, (cf. Kolb & Turner 1990). This result generally obtains after 
recombination and when fluctuations are in the linear regime {dM/M « 1). For example, 
such baryon number perturbations will have grown by a factor of (1 + Zfi) between the 
epoch of recombination and the present epoch if a standard recombination scenario with 
recombination redshift zji ~ 1100 is assumed. 

If the universe stays ionized at redshifts z < zji, the growth of perturbations will be 
inhibited by the coupling of photons to baryons. In the limit where the photon mean free 
path, , is larger than the fluctuation size, I / , growth of baryon number perturbations 
is inhibited by Thomson drag between photons and electrons. This drag force effectively 
suppresses any perturbation growth for redshifts above 200 - 300 (cf. Peebles 

1971). For redshifts below Zdrag, Thomson drag rapidly becomes inefficient and so cannot 
hinder the growth of perturbations. In the limit when < If, sub- horizon scale entropy 
perturbations behave like oscillating sound waves. Note that the photon mean free path 
has a comoving size of roughly ~ 8h~^Mpc (corresponding to Mg) at a redshift of ^ ~ 200. 

In order that cosmic structure in primordial isocurvature baryon number fluctuation 
models not be "overproduced" , it is necessary that the primordial root-mean-square mass 
fluctuations on the mass scale Mg be smaller than {1 + z)~^, 



where 1100^^^200. Here the range in redshifts results from two extreme scenarios: stan- 
dard recombination; and ionization fully maintained down to low redshifts. 

Note that equality in equation (16) pertains to the case when the primordial isocurva- 
ture baryon number fluctuations are the seeds for the presently observed cosmic structure. 
It is, however, also conceivable that the formation of structure on large scales is mainly 
due to, for example, adiabatic fluctuations. In this case perhaps only the small-mass scale 
fluctuations are dominated by the primordial isocurvature baryon number fluctuation com- 
ponent. The inequality in equation (16) would apply in this latter example. 




(16) 
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By using equations (15) and (16) we can derive a limit on the mass cutoff, 

M,<6xl0l3Mo(^^)/^-3{(l + ^^)a,(M,)}"^. (17) 

Our nucleosynthesis results show that this mass cutoff must not fall below Mj. If it did, 
we would produce unacceptable light element abundances. Furthermore, unless ar){Mc) ^ 1 
there would be essentially no interesting effects on nucleosynthesis (e.g., no early collapse 
of overdense regions with significant mass fractions). For a given spectral index it is not 
always clear that these two requirements are not mutually exclusive. 

As an example of a scenario that does meet both criteria consider Mc ~ 3 x 10^ Mq, 
an{Mc) ~ 1, n = 0, and Zr < 250. By contrast, an example of a scenario which will not 
work has ^ 1 x IO^Mq, c7^(Mc) fa 10, n = -1, and Zr ^ 1100. 

2.4. Convergence of numerical results 

We have performed detailed primordial nucleosynthesis calculations with the numer- 
ical techniques outlined above. In Figures 4a and 4b we present several measures of con- 
vergence accuracy as functions of numbers of Fourier modes employed in our numerical 
study. 

For the results shown in Figure 4a we employ < r] >— 6 x 10~^*^, Acr = 1-5, n = 0, 
and the transformation in equation (8a). In Figure (4a) the upper panel shows the ratios 
of ^H/H (solid line), Yp (dotted line), and ^Li/H (short-dashed line) to their converged 
values as functions of the number of zones employed. Also shown in this figure is the ratio 
of Q})/Qdiff to its converged value. Here 0,diff is the fractional contribution of baryons in 
surviving, uncoUapsed, regions to the closure density. It is clear that good convergence is 
obtained for cases with more than a few thousand Fourier modes. 

The lower panel of Figure 4a shows cr^(Mc) (solid line) and 20 x (7^,(1000 x Mc) (dotted 
line) as functions of the number of Fourier modes. Convergence of 0",^ is good whenever 
more than a few thousand Fourier modes are employed. 

For the computations presented in Figure 4b we employ < r] >— 1.2 x 10~^, Acr — 2.0, 
n = 0, and the transformation in equation (8b). The quantities plotted in these figures 
are the same as in Figure 4a, and the notation is the same. Again, we note the rapid 
convergence of light element abundance yields and Q^/Qdiff with increasing number of 
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Fourier modes (upper panel). In the lower panel, however, we must conclude that the 
calculation for cr^ remains unconverged, even for 10^ Fourier modes. 

The reason for the lack of convergence in cr^ for the case in Figure 4b stems from 
the fact that the calculations of cr^ are dominated by (rare) high-cr events in the gaussian 
random variable. Note that this is not the case for the calculations of the light element 
abundances. Therefore, we can predict the light-element abundances in this case with 
confidence, though we are unable to predict accurately the effective spectral index, n , 
which correspond to our distribution for 77 in the cases where the parameters of Figure 4b 
are adopted. Similar conclusions would obtain if we had employed the transformation law 
in equation (8c) instead of that in equation (8b). 

3. Results 

In this section we describe the results of our numerical study of primordial nucleosyn- 
thesis in the presence of large mass scale, PIB-like, entropy fiuctuations. Here we will 
discuss not only the light element (^H, ^He, ^He, ^Li) nucleosynthesis yields from various 
models for the distribution of 77, but also the fraction of baryons in these models which 
survive collapse and comprise the "diffuse" primordial gas. We will characterize the sur- 
viving diffuse component of baryons by its fractional contribution to closure, Qdiffh"^- In 
what follows Qbh"^ refers to the baryonic fraction of the closure density prior to freeze-out 
from nuclear statistical equilibrium and, thus, prior to any significant amount of collapse. 

In Figure 5a and Figure 5b we present nucleosynthesis results for a model of the 
distribution of rj which is characterized by the transformation in equation (8a), spectral 
index n = (employed in equation 7), > Mj , and various values of fift/i^ and Ag^. In 
Figure 5a the panel at upper left gives the ^He mass fraction, 1^, as a function of Q^h^. 
The panel at upper right in this figure gives the ratio of produced to hydrogen, ^H/H, 
as a function of Qbh^. Similarly, ^Li/H and ^He/H versus Qbh^ are shown in the panels 
at lower left and right, respectively. In Figure 5b we plot Qdiffh^ versus Qbh^ for various 
models. In Figure 5a the solid lines give the results for standard HBBN at the indicated 
values of n^/i^. 

The other lines in Figures 5a and 5b are as follows: The dotted line is for Acr = 1-25; 
the short-dashed line is for Acr = 1-5; the long-dashed line is for Acr = 1-75; while the 
dashed-dotted line is for Acr = 2. 
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Note that several general trends are evident in Figures 5a and 5b. First we note that 
for all values of Agr and Q^h'^ considered here Yp is lower than the HBBN yield, while 
^H/H and ^He/H are higher than the HBBN yield. Our models clearly retain the well 
known "^Li dip" , but in general ^Li/H can be slightly smaller, comparable to, or larger 
by a factor up to 5, than the yield from HBBN at a given O^/i^ 

These results are easily understood by comparison of the panels in Figure 5a to the 
graph in Figure 5b. Clearly ^H/H is high in inhomogeneous models compared to HBBN 
at the same Q^h'^ because the surviving diffuse baronic component is characterized by 
^diffh^ < ^bh^- In fact, since in HBBN ^H/H yields rise very steeply with decreasing r] , 
we can identify two competing effects in ^H production. Regions with smaller 77 produce 
relatively more ^H, but they make a relatively smaller contribution to the total surviving 
mass of baryons. We find that the regions which are most effective in producing ^H have 
?7 3 X lO"-'^-'^. Since in primordial nucleosynthesis ^He is produced by ^H(p,7)^He, it is 
not surprising that ^He, like ^H, is always high compared to HBBN at the same Clhh'^. 

The situation for ^He is straightforward. Since Yp is a rising function of 77 in HBBN, it 
is obvious that when Qdiffh'^ < ^bh^ inhomogeneous models will produce low Yp relative 
to HBBN at the same Ob/i^. 

By contrast, the behavior of the ^Li/H yield in our inhomogeneous models is more 
complicated. In part, this is due to the two principle production channels for ^Li. These 
are ^H(Q!,7)^Li, which is dominant in low density regions where ^H (thus also ^H) is high, 
and ^He(Q;,7)^Be(e~,z/e)^Li, which is dominant in higher density regions. Since, to some 
extent, our models average over regions which have baryon-to-photon ratios on opposite 
sides of the dip, our ^Li/H yields are high compared to HBBN for a fair range of Qbh^- 
Note, however, two interesting features of our "^Li/H. results: (1) the dip is offset to higher 
Qbh^ than the HBBN dip; and (2) '^Li/H at higher Qfo/i^ is lower than the yield in HBBN. 

A surprising feature of Figure 5a is that the light element abundance yields in our 
inhomogeneous model are actually fairly insensitive to the values of the parameter Acr 
which is employed. We conclude that, even in cases where both the criterion in Acr 
required for collapse and the collapse efficiencies are not accurately determined, we can 
predict the nucleosynthesis for a given PIB-like fluctuation spectrum with fair confidence. 
This result is quite important since the collapse efficiencies may well depend not only on 
the fluctuation's overdensity but also on internal geometry, intrinsic angular momentum, 
and the local environment (e.g., neighboring fluctuations). We therefore expect collapse 
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efficiencies to not change discontinously from zero to one at some value Ag^, but rather 
approach unity in a continous fashion over some interval of overdensities. 

However, these conclusions are dependent on the assumption that collapsing regions do 
not explode and produce large amounts of explosive nucleosynthesis (c/., Fuller, Woosley, & 
Weaver 1986). Furthermore, the results in Figure 5 depend critically on having a significant 
fraction of the baryons collapse. In the case where only a small fraction of the baryons 
collapse (i.e. Acr very large) we would find Yp and (^Li/H) to be increased considerably 
over the results displayed in Figure 5a. 

At this point we wish to summarize the observational constraints on the light-element 
abundances. The situation for the ^He mass fraction, 1^ , is as follows. The best deter- 
mination of Yp is obtained from observations of ^He-recombination lines in metal-poor, 
extragalactic H H regions. The existing data has been analyzed by a number of authors 
(Walker et al. 1991; Fuller, Boyd, & Kalen 1991; Pagel et al. 1992; Mathews, Boyd, 
& Fuller 1993, Olive & Steigman 1994). It is generally believed that the upper limit on 
Yp should be somewhere in the range Yp < 0.24 — 0.245, with the "favored" value for Yp 
around Yp ^ 0.23. 

However, more recently it has been pointed out that the observational determination 
of Yp is subject to several systematic uncertainties which have previously not been well 
appreciated (Skillman & Kennicut 1993; Skillman, Terlevich, & Garnett 1994; Sasselov & 
Goldwirth 1994). When all the systematic uncertainties are taken into account, such as 
uncertainties associated with the determination of emissivities and corrections due to the 
possible existence of some neutral helium in H II regions, a firm upper limit on Yp may 
be as large as 1^ < 0.252. A lower limit on Yp, which is of less constraining power for 
PIB-like models than the upper limit, should be somewhere around Yp > 0.21 — 0.22. It is 
clear that an accurate determination of the primordial ^He-mass fraction from the existing 
data is actually not that straightforward. 

The situation is not very different for the primordial abundances of ^H, and ^He. It 
is common practice to give the upper limit for the sum of ^H and "^Hc (Walker et al. 
1991). This is done since it is not known to what extent ^H, the most fragile of the light 
elements, has been destroyed in stars prior to the formation of the solar system. Before 
the very recent claims of an extragalactic observation of ^H by Songaila et al. (1994), the 
deuterium-to-hydrogen number ratio had been estimated only for the solar system. By 
considering the sum of ^H and ^He some of the uncertainties in the determination of the 
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individual abundances are evaded, since the destruction of is expected to lead to the 
production of ^He via ^H(p,7)^He. In any case, it is commonly believed that the upper 
limit on the sum of these light elements is C^R+^Re) /R <10-^ (Smith et al. 1993). The 
lower limit for deuterium is usually given at (^H/H) J^l.8 x 10~^. 

Recently, Songaila et al. (1994) observed an isotope-shifted Lyman-cu absorption line 
in a Lyman-a cloud system which lies along the line of sight to a quasar. The existence 
of an absorption line is well explained if the Lyman-o; cloud has a deuterium-to-hydrogen 
number fraction as large as (^H/H) 1.9 x 10~^ — 2.5 x 10~^. This observation may cast 
doubt on the previously believed upper limit on primordial (^H+^He) /H. It may, however, 
also be that the detected absorption line comes from hydrogen as opposed to deuterium. 
This could occur if, by coincidence, a small component of the Lyman-a cloud system has 
a small collective-velocity relative to the main cloud. In such a case, the existence of 
an isotope-shifted Lyman-a absorption line might be mimicked. Further observations of 
Lyman-CK cloud systems will hopefully provide us with a reliable primordial deuterium 
abundance. 

The situation for the primordial ^Li abundance is controversial as well. Spite & Spite 
(1982) detected a ^Li abundance plateau for hot, low-metallicity Population II halo stars 
over some temperature and metallicity range. A commonly held view is that the ^Li- 
abundance of (^Li/H) « 1.4 X 10"-*^° observed in the "plateau" -stars is the primordial 
one (Spite & Spite 1982; Thorburn 1994). When uncertainties in the employed model 
atmospheres for the plateau-stars and "small" amounts of diffusion-induced ^Li depletion 
are taken into account, the upper limit on the primordial value of ^Li may be larger by a 
factor which is roughly between one and two (Deliyannis, Demarque, & Kawaler 1990). 

The above interpretation of the data precludes the possibility of significant ^Li de- 
pletion in plateau-stars. However, it has been shown by Pinsonneault, Deliyannis, & 
Demarque (1992) and Chaboyer & Demarque (1994) that ^Li in Population II stars is 
depleted by up to an order of magnitude, or more, over the entire surface temperature 
range (including the Spite-plateau) when microscopic diffusion and rotation of stars are 
included in the stellar models. Furthermore, only combined models which include rotation 
and diffusion can simultaneously explain the very different ^Li depletion patterns in old 
Population II stars and young Population I stars (Pinsonneault et al. 1992; Chaboyer, 
Demarque, & Pinsonneault 1994ab). This may be a powerful argument for the validity of 
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the combined models. Combined models predict a ^Li-abundance of (^Li/H) ~ 1 x 10~^, 
much larger than the value of the Spite-plateau. 

It has been pointed out that a possible detection of ^Li in two population II stars 
(Smith, Lambert, & Nissen 1982; Hobbs & Thorburn 1994) may provide an argument 
against significant depletion. This is because ^Li, presumbly produced by cosmic ray 
spallation of heavier elements over the entire history of the universe, would have been 
depleted below detection level along with ^Li. On the other hand, the observations of three 
highly ^Li-depleted plateau stars (Thorburn 1994) indicates that significant ^Li-depletion 
in plateau stars may occur, at least for some stars. 

The situation remains controversial. It should be noted that the primordial ^Li abun- 
dance, when determined with confidence, could be a crucial argument against, or indication 
for, the existence of small-scale and/or large-scale baryon-to-photon inhomogeneity during 
the nucleosynthesis epoch. 

The implications of the observatially inferred primordial abundance constraints for the 
model results displayed in Figures 5a and 5b are as follows. There seem to be two possible 
ranges in O^/i^ which could yield agreement with abundance constraints within the above- 
mentioned observationally inferred primordial abundance uncertainties. For O^/i^ ^ 0.015 
wc find from Figure 5a that (^H/H) ^ 2 x 10"^ - 4 x 10"^; Yp ^ 0.225 - 0.235; and 
(^Li/H) 3.5 X 10~^° — 4.5 x 10~^° depending on the employed Acr- To be correct, this 
scenario would have to assume that there was a modest amount of ^Li depletion in halo 
stars and a high primordial deuterium abundance. 

For the range flih^ ^ 0.03 - 0.05, we obtain (^H/H) 4 x 10"^ - 1.5 x 10"^; 
Yp fa 0.24 - 0.25; and (^Li/H) 5 x 10"^° - 1.5 x 10"^ This model would then require a 
significant ^Li depletion, a relatively low deuterium abundance and high ^He abundance. 
For comparison, the preferred range for Q^h^ in homogeneous big bang nucleosynthesis 
is somewhere in the range Qth^ ~ 0.005 — 0.015, mostly depending on what primordial 
deuterium abundance is adopted. 

Another generic feature of our models with large-scale inhomogeneity in i] is evident 
from Figure 5b. Even though the abundance yields in these models tend to agree with 
observationally inferred abundance limits for higher O^/i^ than the allowed range of this 
quantity in standard homogeneous big bang models, the fractional contribution of the 
difi^use, "surviving" baryons to the closure density, Qdiffh^i tends to be lower than the 
allowed range of flbh'^ in homogeneous models. For a good agreement between abundance 
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yields and observationally inferred abundance limits in PIB-like models, the low-density 
regions must have approximately the preferred average standard homogeneous baryon-to- 
photon ratio (i.e. Qbh'^ ^ 0.0125 for (2H/H)<10-^). Therefore, Vtdiffh'^ is lower than 
the Jlft/i^ from a homogeneous big bang, since only the low-density regions in the universe 
contribute to the diffuse baryons. The high-density regions do not contribute since these 
presumbly collapse and form condensed, dark objects. 

This way of lowering Vlcuffh? compared to Vthh? inferred from homogeneous big bang 
nucleosynthesis is quite analogous to the scenario proposed by Jedamzik, Mathews, & 
Fuller (1994). In that paper it was shown that there is no lower limit on Q,i,h?' in inho- 
mogenous primordial nucleosynthesis scenarios, since it may be that only a certain frac- 
tion of space has baryons with baryon-to-photon ratio 77 fa 3 x 10"-^°, with the remaining 
fraction of space depleted in baryons. Such a scenario could provide an explanation for 
the small amount of baryons observed in luminous form (i.e., galaxies and diffuse inter- 
galactic gas). The fractional contribution of luminous baryons to the closure density is 
^lum ^ 0.003 - 0.007. This value can be much smaller than the inferred from homo- 
geneous big bang nucleosynthesis. It is interesting to note that inhomogeneous PIB-like 
models conceivably could provide a natural explanation for the small value of Q}^"^. 

In Figure 6a and 6b we present nucleosynthesis results for a model of the distribution 
of ?7 which is characterized by the transformation in equation (8b), spectral index n = 0, 
Mc> M), and various values of Q.^h'^ and Acr- The notation is as in Figure 5a and 5b. 
The assumed values for A^^ in these calculations are as follows: the dotted line is for 
Ac^ = 1.5; the short-dashed line is for Acr = 2; the long-dashed line is for Ag^ = 2.5; and 
the dashed-dotted line is for Acr = 3. 

It is seen from Figure 6a and 6b that most of the general trends observed in the PIB- 
like model of Figure 5a and 5b are retained. It is evident that there is a range in ilb/i^ 
(0.04^f26/i^^0.07) for which all the abundance constraints may be satisfied. In this range 
we find (^H/H) ^ 9x10-^-3x10"^; Yp ^ 0.235-0.248; and (^Li/H) ^ 7xl0-^°-2xl0-9, 
depending on the value of Acr which is employed. Clearly, such models would be ruled 
out if the ^Li abundance of the Spite-plateau in halo stars represents the actual primordial 
abundance. Low values for r2ft/t^^0.02 are ruled out by deuterium overproduction. While 
the models of Figure 5 are characterized by collapse ratios Qh/^diff ~ 2.5 — 10, a much 
larger fraction of the baryons would collapse in the model of Figure 6, Qb/fldiff ~ 25 — 75. 
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In Figure 7a and 7b we present the results for a distribution in rj caracterized by the 
transformation in equation (8c), spectral index n = and Mg > Mj. As in Figure 5 and 
6 we give results for various fife/i^ and Acr- The lines are for Acr = 1.5 (dotted line), 
Acr = 2 (short-dashed line), Acr = 2.5 (long-dashed line), and Acr = 3 (dashed-dotted 
line) . 

In this model, the low "^He-mass fraction, and high (^H/H)-number ratio compared to a 
HBBN model at the same Qbh^ is most pronounced. The collapse ratios, flb/fldiff^ are very 
similiar to the collapse ratios in the model of Figure 6, in particular flb/Qdiff ~ 25 — 75. In 
the range Qbh^ ~ 0.06 - 0.2 we obtain abundance yields of (^H/H) ft! 4 x 10~^ - 4 x 10~^, 
Yp ft 0.225 - 0.25, and (^Li/H) ft 9 x 10"^° - 3 x 10"^ depending on the value for Acr. 
These abundance yields may agree with observationally inferred abundance limits when 
significant ^Li depletion in population II halo stars occurs. It is remarkable that the allowed 
(albeit with high ^Li) range of flbh^ in these models is between a factor of ten and fourty 
larger than the allowed range of Q^h^ in HBBN models. The upper end of this range in 
Qhh^ may even allow for baryons to provide closure density when the Hubble parameter 
is smaller than h^OAb. 

We have also investigated the nucleosynthesis results in PIB-like models as a function 
of the collapse efficiency parameter / as defined in Section 2. In Figures 8a and 8b we 
present the nucleosynthesis yields and ratios Vtb/^diff for a distribution in ry characterized 
by the transformation in equation (8a), spectral index n = 0, Mc > Mj, and Acr = 1.5. 
In this figure we vary Q,i,h? and the collapse efficiency parameter. The notation in these 
figures is similiar to that in Figures 5,6 and 7. 

The different lines correspond to the following collapse efficiencies: / = 100% (lower 
dotted lines), / = 97.5% (short-dashed lines), / = 95% (long dashed lines), / = 90% 
(dashed-dotted lines) , and / = 80% (upper dotted lines) . Note that the lower dotted lines 
in Figures 8a and 8b represent the results for the same model as the short-dashed lines in 
Figures 5a and 5b. These lines are shown for comparison. Note that the lower dotted lines 
in the Yp and (^Li/H) panels of Figures 8a and 8b are to be associated with the upper 
dotted lines in the (^H/H) and ('^He/H) panels of Figure 8a. 

It is not surprising to find that the ^He and ^Li abundances increase and the ^H and 
^He abundances decrease when the collapse efficiency parameter decreases. This is because 
with lower collapse efficiency more high-density regions contribute their nucleosynthesis 
yields to the diffuse baryons. The high-density regions produce relatively larger amounts 
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of ^He and ^Li and smaller amounts of and ^He than the low-density regions. In this 
model, however, it is evident that collapse efficiencies do not have to be extremely close to 
100% in order to avoid gross overproduction of ^He and ^Li. 

In Figures 9a and 9b we present the nucleosynthesis yields and survival ratios Qi/Qdiff 
for a distribution in i] characterized by the transformation equation (8b), spectral index 
n = 0, Mc > Mj, and A^r = 2. We show models with collapse efficiencies / = 100% 
(dotted lines), / = 99% (short-dashed lines), and / = 97.5% (long-dashed lines). Note 
that the dotted lines in Figures 9a and 9b represent the results for the same model as the 
short-dashed lines in Figure 6a and 6b. 

From Figures 9a and 9b it is evident that the collapse efficiencies for these models 
would have to be very close to 100% in order to avoid overproduction of ^He and, especially, 
^Li. This conclusion stands in contrast to that derived from the models of Figures 8a and 
8b, and is easily understood by an examination of the results presented in Figure lb. This 
figure shows a small part of the distribution in 77 from which the model results of Figures 
9a and 9b have been computed. The distribution is characterized by very overdense peaks. 
If only a very few of these peaks did not collapse, their contributions to the ^Li component 
in the surviving diffuse baryons from these high-density regions would be very significant, 
perhaps even dominant. 

Our computation of the abundance yields from a particular model as a function of 
the collapse efficiency parameter / proceeds in the following manner. First, the algorithm 
finds all the regions which are overdense by the critical amount Acr- Then, the algorithm 
randomly chooses a fraction (1 — /) of these overdense regions to not collapse and thus to 
contribute to the abundance yields of the surviving diffuse baryons. 

This procedure may actually not lead to a very accurate assessment of the variations 
in abundance yields associated with the uncertainties of collapse of particular regions. Re- 
alistically, an overdense region may not completely collapse because of its peculiar angular 
momentum and local environment. However, we expect the very overdense peaks within 
an overdense region to collapse with high efficiency. A better assessment of the variations 
in abundance yields due to collapse efficiency uncertainties may be obtained by averaging 
over the abundance yields for a range of Acr- Figures 5-7, which show the abundance 
yields of different models for different Acr, in^iy therefore provide a more realistic estimate 
for the anticipated magnitude of variations in abundance yields due to collapse efficiency 
uncertainties. In any case. Figure 9a illustrates that a very large fraction (^99%) of the 
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overdense peaks of a distribution characterized by very overdense peaks must collapse in 
order that ^He and ^Li not be overproduced. 

In Figure 10a and 10b we show the nucleosynthesis results and ratios ilf^/fldiff for 
a distribution in rj characterized by transformation equation (8c), spectral index n = 0, 
Mc > Mj, Acr = 2, and varying O^^^ and collapse efficiencies /. The lines represent 
/ = 100% (dotted lines), / = 99% (short-dashed Hues), and / = 97.5% (long-dashed 
lines). The conclusions drawn from these figures are quite similiar to the conclusions for 
the model investigated in Figures 9a and 9b. A large fraction (/^99%) of very overdense 
peaks has to form dark remnants in order that ^Li and ^He not be overproduced. 

In our study we have so far assumed that the fluctuation cutoff mass scale, Mg, is 
larger than the mass scale Mj. The local baryon Jeans mass, Mj, divides fluctuation 
evolution into two regimes: overdense fluctuations on mass scales M > Mj ultimately are 
expected to collapse; whereas, overdense fluctuations on mass scales M < Mj will expand, 
and their nucleosynthesis yields will mix with the diffuse baryons. We have examined the 
nucleosynthesis yields of PIB-like models when fluctuations exist on mass scales below Mj. 

In Figures 11a and lib we show the nucleosynthesis yields and fraction Qh/Qdiff for 
a PIB-like model with a distribution in r] characterized by the transformation equation 
(8a), spectral index n — 0, Acr = 1-5, and collapse efficiency / = 100%. In this figure we 
have varied the cutoff mass scale Mc- The dotted line is for Mc > Mj, the short-dashed 
line is for Mc = Mj/3, the long-dashed line is for Mc = Mj/6, and the dashed-dotted line 
is for Mc = Mj/ 10. Note that the dotted lines represent the results of our "reference" 
model, which have already been shown by the short-dashed lines in Figure 5a and 5b and 
the dotted lines in Figure 8a and 8b. 

Compared to the reference model, the ^He and ^Li yields increase and the and ^He 
yields decrease whenever Mc falls below Mj. These results are as expected, since sub- Jeans 
mass size overdense fluctuations which are not included in super-Jeans mass size overdense 
regions will expand and contribute large amounts of ^He and ^Li to the diffuse baryons. 
This material from sub- Jeans mass size overdense fluctuations will also dilute the (^H/H) 
and ("^He/H) ratios relative to the reference model. Note that ^He and ^Li (for most of the 
"interesting" range in flbh'^) will even increase over the HBBN yields as the fluctuation 
cutoff mass scale is decreased. 

It is evident from Figure lib that the fraction of material included in regions which 
meet the requirement to have overdensity Acr on a mass scale larger than Mj decreases 
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as the cutoff mass scale decreases. Equivalently, the ratio Q,b/^diff decreases with 
decreasing Mc as well. In the limit where Qb/Qdiff approaches unity, there is no collapse 
and the abundance yields would be given by a weighted average over the HBBN yields 
of all the regions of the entire distribution. The model results displayed by the dashed- 
dotted lines are actually not far removed from that limit {Qi/fldiff^'^)- For comparison, 
the same model but without any collapse {Qh/Qdiff = 1) would yield (^H/H) = 1.5 x 10~^, 
Yp = 0.242, and (^Li/H) = 8.4 x 10"^° for Q^h'^ = 0.01. These abundance yields differ 
only slightly from the yields for Q.^h'^ = 0.01 (dashed-dotted line). 

It is quite surprising that even for a distribution of 77 as inhomogeneous as that in 
Figure la, and without any collapse, there exists a narrow interval in Q^h^ where all the ob- 
servational abundance constraints may be satisfied. In particular, for 0.006;5^6^^^0.0012 
agreement between the computed abundance yields and the observationally inferred abun- 
dance constraints is possible within the observational uncertainties. Note that this range 
for Qijh'^ is close to the range for O^/i^ inferred from HBBN. Agreement between abun- 
dance yields and observationally inferred abundance limits would require significant ^Li 
depletion and a high deuterium abundance (^H/H)^l x 10~^. 

We have implicitly assumed so far that there is a cutoff mass scale which is larger 
than Mc^lO~^^ Mq (approximately sixteen orders of magnitude smaller than Mj\). If 
there were fiuctuations on mass scales below M^10~^^Mq, diffusive and hydrodynamic 
processes during the epoch of primordial nucleosynthesis would alter the abundance yields 
from what we have calculated. 

It is not straightforward to estimate accurately the effects of diffusive processes during 
the nucleosynthesis era on the abundance yields of stochastic baryon number distributions 
such as those shown in Figures labc. The nucleosynthesis yields of a regular lattice of 
fluctuation sites in the mass range 10~'^^ Mq^M'^10~^^ Mq have been investigated in detail 
(Kurki-Suonio et al. 1988; 1990; Mathews et al. 1990; Jedamzik et al. 1994). Depending 
on the fluctuation characteristics, especially the separation of adjacent fluctuations in the 
regular lattice, Yp is usually larger (but can be smaller) than the Yp of a homogeneous model 
at the same average O^/i^. The number ratio of (^Li/H) tends to be much larger in such 
inhomogenous models than its HBBN value. On the other hand, most inhomogeneous 
scenarios with diffusion tend to yield slightly lower values for Yp and (^Li/H) than do 
inhomogeneous models which neglect diffusive processes. 
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The abundance yields for a truly stochastic distribution may be roughly approximated 
by averaging over the abundance yields from different regular lattices of fluctuation sites 
with varying fluctuation separation lengths. Such an average should always increase Yp 
and (^Li/H) relative to their respective HBBN yields (Meyer et al. 1991). Therefore, 
we expect models which have small-scale fluctuations and diffusive processes during the 
nucleosynthesis era to produce ^He and ^Li yields which are far above the dotted lines 
shown in Figure 11a. We also expect such models to yield ^Li and ^He slightly below 
the dashed-dotted lines in Figure 11a. The situation for the other light elements, and 
^He, is more complicated. These issues are investigated in more detail by Kurki-Suonio, 
Jedamzik, & Mathews (1994). These authors will treat diffusive processes during the 
nucleosynthesis epoch explicitly and combine their results with the results of the present 
study. 

Note that our arguments implicitly assume that there is some turnover in effective 
spectral index of the baryon number distribution, from n^ff = on large mass scales, 
to an Ueff = —3 Harrison-Zeldovich character on small mass scales. If there were not 
such a turnover in effective spectral index, there would either be no significant fraction of 
baryons collapsing (since fiuctuation amplitudes on Mj-size mass scales are very small), 
or fluctuation amplitudes on small mass scales M^10~^^Mq would have to be extremely 
large (A>10S). 

In Figures 12a and 12b we show the nucleosynthesis yields and fraction fli,/Qdiff 
for a distribution of rj characterized by the transformation equation (8b), spectral index 
n = 0, Acr = 2, and varying fl^h^ and Mg. The lines represent results for > Mj 
(dotted), Mc = M5/3 (short-dashed), = (long-dashed), and = M^/IO 

(dashed-dotted). In these flgures the dotted line represents the results of our reference 
model which have already been shown by the short-dashed lines in Figures 6a and 6b and 
the dotted lines in Figures 9a and 9b. 

When the cutoff mass scale is decreased, Yp and (^Li/H) abundance significantly 
increase, and (^H/H) and ("^He/H) decrease relative to the results of the reference model. 
Since this particular distribution in ry includes very overdense peaks (A ~ 100), it is not 
surprising that the ^He and ^Li yields increase over the yields of a HBBN scenario at the 
same ri(,/i^. 

For the particular distribution investigated here, there seems to be no interval in 
Qbh^ for which all the abundance constraints may be met. This result implies that models 
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with stochastic large- and small-scale, large-amplitude (A ~ 100) inhomogeneities in the 
baryon-to-photon ratio are ruled out by the observationally inferred primordial light ele- 
ment abundances. We can only expect possible agreement between the abundance yields of 
such models and observationally inferred abundance limits when there is a fluctuation cut- 
off mass scale Mc^Mj. Furthermore, our discussion about the effects of diffusive processes 
during the nucleosynthesis era should indicate that the existence of small-scale fluctuations 
down to the mass range M ~ 10~^^ — 10~^^Mq is not likely to change these conclusions. 

In Figures 13a and 13b we present the nucleosynthesis yields from a distribution in rj 
characterized by transformation equation (8c), spectral index n = 0, Acr — 2, and for cutoff 
mass scales Mc > (dotted), = (short-dashed), = M]/Q (long-dashed), 

and Mc = Mj/ 10 (dashed-dotted). It is evident from these flgures that for < Mj there 
seems to be no range in Clbh'^ where the production of all the light elements is consistent 
with observationally inferred abundance constraints. Note that this model is characterized 
by large-amplitude fluctuations, and so is similiar to the model considered in Figures 12a 
and 12b. 

We have investigated the dependence of nucleosynthesis yields in PIB-like models on 
the effective spectral index of the distribution in r]. We have employed spectral indices 
n = and n = —2.4 in equation (7). Note that n here denotes the three-dimensional 
spectral index. The ar^{M) (or, equivalently, 5p/ p as a function of mass scale) for the 
resultant distributions in -q has been presented in Figures 3abc. These flgures illustrate that 
the effective spectral index of the distribution in -q is rieff ~ for n = and ng// ~ — 1.5 
for n = —2.4. We note, however, that there is some uncertainty associated with the exact 
value of TT-e//. 

In Figures 14a and 14b we show the nucleosynthesis yields and flb/^diff for a dis- 
tribution in rj characterized by transformation equation (8a), Acr = 1.5, Mc > Mj, and 
effective spectral indices rie// ~ (dotted lines) and n^ff ~ —1.5 (dashed lines). It is 
surprising to flnd that there is so little dependence of the nucleosynthesis yields on the ef- 
fective spectral index of the distribution in i]. In general, for decreasing Ueff the ^He mass 
fraction, Yp, decreases slightly, whereas the number ratios (^H/H) and (^He/H) increase 
slightly. The number ratio (^Li/H) decreases for large ri^/t^, and increases for small $1?,/*^ 
when Ucff decreases. 

These trends can be easily understood on inspection of Figure 14b. For a smaller ng//, 
a larger fraction of the baryons collapse. In turn, a larger fraction of collapsing baryons 

27 



implies that the average Q.^iffh? is smaller, which ultimately leads to smaller ^He yields 
and larger and ^He yields. The ^Li yield increases or decreases for decreasing 
depending on whether Q.^h? is smaller than, or larger than, the Q.hh?' at the ^Li dip. Our 
results show that the nucleosynthesis yields in PIB-like models are not very dependent on 
the effective spectral index of the distribution of rj in the interval — 1.5^?ie//^0. Note that 
PIB-models for large scale structure formation would prefer spectral indices somewhere 
between n^ff = and Tie// = —1. 

These conclusions are confirmed by the results shown in Figures 15a and 15b and 
Figures 16a and 16b. These figures show nucleosynthesis yields and Q,i,/Q,diff for models 
with = 2, Mc > Mj, and two different spectral indices: n^ff ~ (dotted lines) 
and rieff ~ —1.5 (dashed lines). The calculations described by Figures 15a and 15b 
employ transformation equation (8b) for the generation of the distribution of 77; whereas, 
calculations described by Figures 16a and 16b employ transformation equation (8c) for the 
generation of the distribution of r}. 

There is a unique prediction of cosmological models which contain non- linear, interme- 
diate to large scale primordial isocurvature baryon number fluctuations. Since such models 
are characterized by inhomogeneity in the baryon-to-photon ratio during the nucleosynthe- 
sis epoch, we expect the production of intrinsic spatial variations in the primordial light 
element abundances. In contrast, HBBN models predict a universal set of cosmic light 
element abundances. 

To illustrate this point, we define the probability distribution function P\{Yp). This 
function gives the probability for finding a region of size A to have average ^He mass 
fraction Yp. Note that the definition of P\iYp) is analogous to the definition of P\{ri) in 
Section 2. In a similar fashion we define probability functions Pa(^H/H), Px(^He/H), and 
P\(^Li/H) to give the probability for finding a region of size A with average number ratios 
of (^H/H), (^He/H), and (^Li/H), respectively. 

In Figures 17a, 17b, and 17c we display these numerically determined probability 
functions for various scales A and different models of the distribution of rj. In these figures 
the panels in the upper left-hand corners display P\{Yp), whereas the panels in the lower 
left-hand corners display P\(JlA/Yi). The panels in the upper and lower right-hand corners 
display Px(^H/H) and Px(^He/H), respectively. In each panel we show four probability 
distribution functions. These distribution functions are determined for the scales A < 
Ac (solid lines), A = 100 Ac (dotted lines), A = 1000 Ac (short-dashed lines), and A = 



28 



lOOOOAc (long-dashed lines). Note that in our "one-dimensional" theory the length scale 
A is proportional to the mass scale M. 

In Figure 17a we show probability distribution functions for a model spatial distribu- 
tion of 1] which is characterized by the transformation equation (8a), spectral index n = 0, 
Acr = 1.5, and Mc > Mj. Furthermore, we assume a cosmic average baryon-to-photon 
ratio of < ry >= 6 x 10" This corresponds to a value of Qbh'^ = 0.0224. 

The panels in Figure 17a clearly illustrate that there is a finite width to the abun- 
dance probability distribution functions. In other words, there is a finite probability to 
find regions of mass scale M to have abundances which are smaller, or larger, than av- 
erage cosmic abundances. The widths of the distribution functions decrease as the mass 
scale (or, equivalently. A) increases. On the smallest mass scales (solid lines), there is a 
fairly wide range of about equally probable primordial abundances. On the largest mass 
scales (long-dashed lines), the intrinsic widths of the probability distribution functions 
are approximately 4 x 10~^ < Yp > for the ^He mass fraction, 10~^ < (^H/H) > for 
the deuterium-to-hydrogen number ratio, 5 x 10"^ < (^He/H) > for the ^He-to-hydrogen 
number ratio, and 7 x 10~^ < (^Li/H) > for the ^Li-to-hydrogen number ratio. Here the 
brackets denote the cosmic average abundances for this particular model for the distribu- 
tion of 77. 

In Figure 17b we show probability distribution functions for a spatial distribution 
of r] which is characterized by the transformation equation (8b), spectral index n = 0, 
Acr = 2, and Mc > Mj. In this model we assume a cosmic average baryon-to-photon 
ratio oi < r] >= 1.2 x 10~^. We observe the same trends and features in this figure as 
were observed in Figure 17a. However, the intrinsic widths in the distribution functions 
for the largest scales (A = 10^ Ac) are much larger than those in Figure 17a. This is 
because the distribution investigated in Figure 17b includes large-amplitude fluctuations 
in ?7, whereas the fluctuation amplitudes of the distribution investigated in Figure 17a are 
moderate. On the largest scales we find approximate widths of the distribution functions 
to be 2 X 10-2 < Yp > for the ^He mass fraction, 0.4 < (^H/H) > for the deuterium- 
to-hydrogen number ratio, 0.25 < (^He/H) > for the ^He-to-hydrogen number ratio, and 
0.3 < (^Li/H) > for the ^Li-to-hydrogen number ratio. 

In Figure 17c we display distribution functions for a spatial distribution of rj which 
is characterized by the transformation equation (8c), spectral index n = 0, Acr = 2, and 
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Mc > Mj. For this calculation we have assumed a cosmic average baryon-to-photon ratio 
of < 7/ >= 3 X 10-9. 

For the smallest scales, the probability distribution functions display a peak which is 
not centered at the cosmic average abundances. This is easily understood by inspection of 
Figure Ic, where it is obvious that in the present model there exists a minimum baryon- 
to-photon ratio rjmin- A large fraction of the universe in this model has rjmim so that the 
probability to find a small region with abundance yields pertaining to ?7mm is large. The 
intrinsic wdths of the distribution functions on the largest scales are somewhat smaller 
than those found in Figure 17b. 

These results imply that PIB-like models which have an intermediate-scale, non-linear 
fluctuation component would lead to intrinsic primordial abundance variations on small, 
as well, as large mass scales. Note that a mass scale of M = lO^Mg assumed for the 
probability distribution functions shown in Figures 17abc can easily correspond to mass 
scales as large as M ~ 10^^ — lO^^M©, depending on the assumed value of the cutoff 
mass scale, Mg. We do not expect small-scale primordial abundance variations to survive 
to the present epoch. This is because mixing mechanisms, such as shock waves induced 
by supernovae explosions, should be efficient enough to mix the primordial material on 
intermediate mass scales. 

It is, however, questionable if mixing could erase preexisting primordial abundance 
variations on mass scales as large as M ~ 10^ — 1O^^M0. In fact, it is well known that the 
^He mass fraction observed in metal-poor, extragalctic H II regions exhibits a significant 
variation between objects of the same metallicity. Typical '*He mass fractions are in the 
range 0.22 — 0.26. It has been suggested that this spread is not due to observational 
uncertainties, but rather represents a real physical variations. However, even if these 
suggestions are correct, there are several chemical evolution effects which may produce 
such a spread (Campbell 1992). 

In any case, the intrinsic spread in the ^He mass fractions observed in metal-poor H 
II regions may be used to put upper limits on the magnitude of preexisting primordial 
abundance variations. In order to derive useful upper limits on large-scale inhomogeneity, 
the details of the mixing of primordial material would have to be considered. In this way, 
cosmological models which include primordial, non-linear, isocurvature baryon number 
fluctuations could be further constrained. 
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4. Conclusions 



We have examined the primordial nucleosynthesis process in the presence of large mass 
scale, nonlinear entropy fluctuations. A variety of fluctuation spectra were considered, 
including some which have the characteristics of the spectra of PIB models extrapolated 
to smaller mass scales. Our computations provide for the collapse of overdense regions 
with masses above the local baryon Jeans mass, Mj (equation 4). The baryons, and hence 
the nucleosynthesis products, which avoid incorporation into condensed objects were found 
to originate mostly from underdense regions. A complicating feature of nucleosynthesis 
calculations in a field of stochastically distributed fluctuations arises from the fact that a 
particular underdense region may reside inside a larger overdense region which, in turn, 
might be destined for collapse. In our computations we have included a detailed numerical 
treatment of this "cloud-in-cloud" problem, and we have found that such a treatment was 
important for the accurate estimate of light element (^H, "^He, ^He, '^Li) abundance yields. 

In general, we have found that PIB-like spectra with significant small-scale, large- 
amplitude structure (i.e., non-linear, A >> 1, structure on mass scales M < Mj) produced 
light element abundance yields which were in confiict with observationally-inferred limits 
for almost any pre-nucleosynthesis Qb- This may represent an important nucleosynthesis- 
based constraint on PIB models for large scale structure formation, if the underlying 
microscopic mechanisms which generated fiuctuations in these models somehow demanded 
the presence of nonlinear small-scale structure. 

However, if such small scale fiuctuations were absent or suppressed, then our compu- 
tations have shown that there exists a range of fiuctuation spectral characteristics which 
produce light element abundance yields in agreement with observationally-inferred lim- 
its. Particular distributions of baryons were found to produce acceptable nucleosynthesis 
yields even when the pre-nucleosynthesis baryonic fraction of the closure density was as 
high as rtb ~ 0.2/1-2 {i.e., roughly closure density for small Hubble parameter). On the 
other hand, the fractional contribution of the diffuse, "uncoUapsed" baryons to the closure 
density, Vtdiff^ was found to tend to be lower than the O™^'^ inferred from standard 
homogeneous big bang nucleosynthesis. Typical values for ^diff < 0.01 may be in better 
accord with the observed fractional contribution of the luminous baryons to the closure 
density {nium ~ 0.003 - 0.007) than VL^^^^ . 

In any case, we have found that such a relaxation of the homogeneous big bang limit 
on fifo would usually require that the primordial abundance of ^Li be closer to the observed 
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Population I value than to the Spite "plateau" value. In turn, this would demand that 
there had been a fair amount of destruction of ^Li in the Spite plateau stars, a conceivable 
though controversial possibility. Future observations may resolve this issue. 

In our models which met abundance constraints we found that ^H/H was high and 
the ^He mass fraction low relative to a homogeneous big bang at a given value of Qbh^- 
Obviously, this would have to be the case if these models were to meet abundance limits 
for values of that were larger than the limit on this quantity from the homogeneous big 
bang. The important point is that PIB-like models possessing small scale fluctuation cutoffs 
could alter the relative abundances of ^H, "^He, ^He, ^Li over those from a homogeneous big 
bang, while still meeting abundance constraints within observational uncertainties. More 
accurate observational determinations of any two of these light element abundances might 
provide a signature for or, more likely, a constraint on such PIB-like models. 

We have pointed out that a potential signature of a PIB-like distribution of entropy at 
the nucleosynthesis epoch would be the observation of a significant and intrinsic variation 
in the primordial abundances of the light elements, especially deuterium. Our calculations 
have predicted that such PIB-like models could produce light element abundance variations 
on baryon mass scales up to 10"'^° M© to lO"*^^ M0. 
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6. Figure Captions 

Figure la A one-dimensional baryon-to-photon number distribution r]{x) divided by the average 
baryon-to-photon ratio < ry > as a function of space coordinate x (solid line). This 
distribution has been generated by employing the transformation equation (8a) to a 
gaussian random variable. The dotted boxes indicate those regions which are over- 
dense on average by the critical amount Acr- In this figure we have used Acr = 1-5. 
Note that the full simulation- "volume" extends from a; = to a; = 1. The figure shows 
only a small part of the full distribution. 

Figure lb The notation in this figure is as in Figure la. The distribution shown in this figure has 
been generated by employing the transformation equation (8b) to a gaussian random 
variable. For this figure we have used a critical overdensity of Acr = 2. 

Figure Ic The notation in this figure is as in Figure la. The distribution shown in this figure has 
been generated by employing the transformation equation (8c) to a gaussian random 
variable. For this figure we have used a critical overdensity of Acr = 2. 

Figure 2a Probability distribution functions P\{r]) to find a region of size A (or, equivalently, 
mass M) with average baryon-to-photon ratio 77. These distribution functions are 
plotted as a function of the ratio r]/ < rj >, where < r] > is the cosmic average baryon- 
to-photon ratio. Distribution functions are shown for the mass scales M = (solid 
line), M = lOMc (dotted line), M = lOOMc (short-dashed line), and M = lOOOMc 
(long-dashed line). Here is the cutoff mass scale below which baryon-to-photon 
fiuctuations are assumed to be suppressed. For comparison, the dashed-dotted line 
shows the analytically determined distribution function P{r]) (equation 11). Note 
that the ordinate scales differently for different distribution functions. The baryon- 
to-photon distribution has been generated by employing the transformation equation 
(8a) to a gaussian random variable. 

Figure 2b The notation in this figure is as in Figure 2a. Distribution functions are shown for 
a baryon-to-phton distribution generated by employing transformation equation (8b) 
to a gaussian random variable. 

Figure 2c The notation in this figure is as in Figure 2a. Distribution functions are shown for a 
baryon-to-photon distribution generated by employing transformation equation (8c) 
to a gaussian random variable. 

Figure 3a The variance cr^(iW^) as a function of the ratio of mass to cutoff mass scale [M/Mc). 

We have generated the baryon-to-photon distribution by employing transformation 



36 



equation (8a) to a gaussian random variable. The dotted line shows cr^(M) for a 
distribution where a spectral index n = of the gaussian random variable has been 
used (equation 7). The dashed line shows (T^(M) for a model where n = —2.4 has 
been used. For comparison we also show (jn{M) from equation (15) with n = (lower 
solid line) and n = —1.5 (upper solid line). 

Figure 3b The notation in this figure is as in Figure la. For this figure we use a distribution in 
T] which is characterized by transformation equation (8b). 

Figure 3c The notation in this figure is as in Figure 3a. For this figure we use a distribution in 
T] which is characterized by transformation equation (8c). 

Figure 4a Convergence of numerical results as a function of the number of Fourier modes em- 
ployed in the simulation. Upper panel - the ratios (^H/H) (solid line), '^He mass 
fraction Yp (dotted line), (^Li/H) (short-dashed line), and (Qi,/Qdiff) to their con- 
vergent values as noted in the figures. Lower panel - the variances a"^(Mc) (solid 
line), and 20 x (j^(1000 x M^) (dotted line). In these simulations we have used the 
transformation equation (8a), spectral index n = 0, and Acr = 1-5. 

Figure 4b The notation in this figure is as in Figure 4a. For this figure we have used a distribution 
in T] characterized by transformation equation (8b), spectral index n = 0, and Acr = 2. 

Figure 5a Light-element nucleosynthesis yields for a model distribution of r] which is charac- 
terized by the transformation equation (8a), spectral index n = 0, Mc > Mj, and 
various values of O^^^ and A^- The panel in the upper left-hand corner shows the 
"^He mass fraction 1^ as a function of Vtbh? ^ whereas the panel in the lower left-hand 
corner shows the (^Li/H) number ratio as a function of Vti,h? . The panels in the up- 
per and lower right-hand corners show the (^H/H) and (^He/H) number ratios as a 
function of ilfo/i^, respectively. The dotted line is for Acr = 1-25, the short-dashed 
line is for Acr = 1-5, the long-dashed line is for Acr = 1-75, while the dashed-dotted 
line is for Acr = 2. The solid lines give the results of standard homogeneous big bang 
nucleosynthesis for comparison. 

Figure 5b The values for Q^iffh'^ as a function of fli,h'^ for those models for which the nucle- 
osynthesis yields are shown in Figure 5a. 

Figure 6a Nucleosynthesis yields for a model distribution of rj which is characterized by the 
transformation equation (8b), spectral index n = 0, Mc > Mj, and various flbh'^ and 
Acr- The different lines show results for models with Acr = 1-5 (dotted line), Acr = 2 
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(short-dashed hne), Agr = 2.5 (long-dashed hne), and Agr = 3 (dashed-dotted hne). 
The notation in this figure is as in Figure 5a. 

Figure 6b The values for Q.^iffh? as a function of Q.hh? for those models for which the nucle- 
osynthesis yields are shown in Figure 6a. 

Figure 7a Nucleosynthesis yields for a model distribution of rj which is characterized by the 
transformation equation (8c), spectral index n = 0, Mc > Mj, and various O^/i^ and 
Acr- The different lines show results for models with Acr = 1.5 (dotted line), A^r = 2 
(short-dashed line), Acr = 2.5 (long-dashed line), and Acr = 3 (dashed-dotted line). 
The notation in this figure is as in Figure 5a. 

Figure 7b The values for Qdiffh'^ as a function of Q^h'^ for those models for which the nucle- 
osynthesis yields are shown in Figure 7a. 

Figure 8a Nucleosynthesis yields for a model distribution of r] which is characterized by the 
transformation equation (8a), spectral index n = 0, Mc > Mj, Acr = 1-5, and various 
fifo/i^ and collapse efficiency parameters /. The different lines show results for models 
with / = 100% (lower dotted line in the panels for Yp and (^Li/H)), / = 97.5% 
(short-dashed line), / = 95% (long-dashed line), / = 90% (dashed-dotted line), and 
/ = 80% (upper dotted line in the panels for Yp and (^Li/H)). Note that the lower 
dotted lines in the panels for Yp and (^Li/H) correspond to the upper dotted lines in 
the panels for h2 and ("^He/H), and vice versa. The notation in this figure is as in 
Figure 5a. 

Figure 8b The values for Q^iffh^ as a function of Q^h^ for those models for which the nucle- 
osynthesis yields are shown in Figure 8a. In this figure the results shown by the lower 
dotted line correspond to the results shown by the lower dotted lines in the panels for 
Yp and C^Li/H) in Figure 8a. 

Figure 9a Nucleosynthesis yields for a model distribution of r] which is characterized by the 
transformation equation (8b), spectral index n = 0, > Mj, Acr = 2, and various 
flbh'^ and and collapse efficiency parameters /. The different lines show results for 
models with / = 100% (dotted line), / = 99% (short-dashed fine), and / = 97.5% 
(long-dashed line). The notation in this figure is as in Figure 5a. 

Figure 9b The values for VtdiffhF' as a function of VL^fh? for those models for which the nucle- 
osynthesis yields are shown in Figure 9a. 
Figure 10a Nucleosynthesis yields for a model distribution of 77 which is characterized by the 
transformation equation (8c), spectral index n = 0, Mc > -^j, Acr = 2, and various 
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Qhh^ and collapse efficiency parameters /. The different lines show results for models 
with / = 100% (dotted line), / = 99% (short-dashed line), and / = 97.5% (long- 
dashed line). The notation in this figure is as in Figure 5a. 

Figure 10b The values for Qdiffh'^ as a function of Qbh'^ for those models for which the nucle- 
osynthesis yields are shown in Figure 10a. 

Figure 11a Nucleosynthesis yields for a model distribution of i] which is characterized by the 
transformation equation (8a), spectral index n — 0, Acr = 1-5, and various O^/i^ and 
cutoff mass scales Mc- The different lines show results for models with Mc > Mj 
(dotted line), = (short-dashed line), Mc = (long -dashed line), and 

Mc = Mj/10 (dashed-dotted line). The notation in this figure is as in Figure 5a. 

Figure lib The values for Cldiffh'^ as a function of ilfe/i^ for those models for which the nucle- 
osynthesis yields are shown in Figure 11a. 

Figure 12a Nucleosynthesis yields for a model distribution of r] which is characterized by the 
transformation equation (8b), spectral index n = 0, Acr = 2, and various Qi,h^ and 
cutoff mass scales Mc. The different lines show results for models with Mc > Mj 
(dotted line), Mc = M^/S (short-dashed line), Mc = M]/Q (long-dashed line), and 
Mc = Mj/10 (dashed-dotted line). The notation in this figure is as in Figure 5a. 

Figure 12b The values for Vtdiffh'^ as a function of VLi,h^ for those models for which the nucle- 
osynthesis yields are shown in Figure 12a. 

Figure 13a Nucleosynthesis yields for a model distribution of r] which is characterized by the 
transformation equation (8c), spectral index n = 0, Acr = 2, and various Vthh"^ and 
cutoff mass scales Mc- The different lines show results for models with Mc > Mj 
(dotted line), Mc = M^/3 (short-dashed line), Mc = M^/6 (long-dashed line), and 
Mc = Mj/10 (dashed-dotted line). The notation in this figure is as in Figure 5a. 

Figure 13b The values for Qdiffh^ as a function of flbh^ for those models for which the nucle- 
osynthesis yields are shown in Figure 13a. 

Figure 14a Nucleosynthesis yields for a model distribution of i] which is characterized by the 
transformation equation (8a), Mc > Mj, Acr = 1-5, various O^/i^, and two different 
effective spectral indices Ucf / for the distribution of i]. The dotted line is for a model 
with Ueff ~ 0, whereas the dashed line is for a model with Ueff ~ —1.5. The notation 
in this figure is as in Figure 5a. 

Figure 14b The values for O.diffh'^ as a function of flth'^ for those models for which the nucle- 
osynthesis yields are shown in Figure 14a. 
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Figure 15a Nucleosynthesis yields for a model distribution of rj which is characterized by the 
transformation equation (8b), Mg > Mj, Agr = 2, various Q^h'^, and two different 
effective spectral indices Uef / for the distribution of rj. The dotted line is for a model 
with Ueff ~ 0, whereas the dashed line is for a model with ng// ~ —1.5. The notation 
in this figure is as in Figure 5a. 

Figure 15b The values for Qdiffh'^ as a function of flbh'^ for those models for which the nucle- 
osynthesis yields are shown in Figure 15a. 

Figure 16a Nucleosynthesis yields for a model distribution of i] which is characterized by the 
transformation equation (8c), > Mj, Acr = 2, various Q^h'^, two different effective 
spectral indices Uef / for the distribution of rj. The notation in this figure is as in Figure 
5a. 

Figure 16b The values for fldiffh"^ as a function of ri(,/t^ for those models for which the nucle- 
osynthesis yields are shown in Figure 16a. 

Figure 17a The probability distribution functions Pa(^H/H), PA(^He/H), Px{Yp), and PA(^Li/H) 
as defined in the text. The panel in the upper left-hand corner shows P\{Yp), while 
the panel in the lower left-hand corner shows Px(^Li/H). The panels in the upper and 
lower right-hand corners show Pa(^H/H) and PA(^He/H), respectively. Each panel 
shows four probability distribution functions determined on different scales A. The 
solid line is for A < Ac, the dotted line is for A = 100 Ac, the short-dashed line is 
for A = 1000 Ac, whereas the long-dashed line is for A = 10000 Ac. The distribution 
functions in this figure have been computed from a model spatial distribution in t] 
characterized by the transformation equation (8a), spectral index n = 0, > Mj, 
Acr = 1-5, and average baryon-to-photon ratio < r] >= 6 x lO"^'^. 

Figure 17b The notation in this figure is as in Figure 17a. The probability distribution functions 
shown here have been computed from a model spatial distribution of i] which is char- 
acterized by transformation equation (8b), spectral index n — 0, M^. > Mj, Acr — 2, 
and <ri>= 1.2 x 10"^. 

Figure 17c The notation in this figure is as in Figure 17a. The probability distribution functions 
shown here have been computed from a model spatial distribution of rj which is char- 
acterized by transformation equation (8c), spectral index n = 0, Mc > Mj, Acr = 2, 
and < 7/ >= 3 X 10"^. 
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